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1.  Introduction 


The  split-Hopkinson  pressure  bar  (SHPB)  experimental  technique  has  become  a  standard  method 
in  determining  the  dynamic  mechanical  properties  of  materials.  The  Hopkinson  pressure  bar 
technique  was  originally  developed  by  B.  Hopkinson  (1)  in  1914  to  measure  the  pressure 
generated  by  the  detonation  of  explosives.  In  1948,  R.  M.  Davies  (2)  conducted  a  critical  study 
of  the  Hopkinson  pressure  bar  technique.  He  developed  the  axial  and  radial  strain  measurement 
system  and  analyzed  the  wave  dispersion  effects.  In  1949,  H.  Kolsky  (5)  invented  the  SHPB 
technique.  In  his  work,  he  used  the  same  system  as  Davies’  to  measure  the  strain.  At  the  same 
time,  he  developed  a  one-dimensional  (1-D)  model  to  perform  experimental  data  analysis.  The 
work  of  Hopkinson,  Davies,  and  Kolsky  established  the  foundation  for  the  development  of  the 
SHPB  technique.  Figure  1  shows  a  classic  SHPB  test  setup,  which  consists  of  one  striker  bar 
(SB),  one  incident  bar  (IB)  and  one  transmission  bar  (TB),  equal  in  diameter  and  usually  made 
from  the  same  material.  Strain  gages  are  mounted  on  the  surfaces  of  IB  (SGI)  and  TB  (SG2)  to 
record  the  strain.  The  specimen  (SP)  is  sandwiched  between  the  IB  and  the  TB.  The  SB, 
powered  by  pressured  gas,  will  hit  the  IB.  A  compressive  stress  wave  will  travel  along  the  IB 
and  reach  the  IB-specimen  (IB-S)  interface.  At  this  interface,  one  part  of  it  will  be  reflected  back 
into  IB,  and  the  rest  will  be  transmitted  to  the  TB  through  the  specimen-transmission  bar  (S-TB) 
interface.  The  strain  gage  SGI  will  record  two  waves:  incident  wave  and  reflected  wave  and  the 
strain  gage  SG2  will  record  only  the  transmitted  wave.  Data  reduction  on  these  three  waves  will 
provide  the  stress-strain  response  of  the  specimen  material  under  dynamic  loading. 


Figure  1.  Schematic  diagram  of  SFtPB  setup. 


The  classical  SHPB  data  analysis  is  based  on  1-D  wave  propagation  theory  in  long  rods.  Several 
assumptions  have  been  made  for  the  1-D  theory.  Two  important  assumptions  are  as  follows: 

(1)  stress  wave  propagation  in  the  bars  is  1-D,  and  (2)  the  bars  are  free  of  dispersion  effects. 
Actually,  these  two  assumptions  are  not  satisfied  for  finite  diameter  bars.  Based  on  such 
assumptions,  the  strain  data  measured  at  the  middle  point  of  the  bars  are  directly  used  as  the 
strain  at  the  interfaces  to  reduce  the  data.  Strictly  speaking,  wave  propagation  in  cylindrical  bars 
is  three-dimensional  (3-D)  in  nature:  both  axial  and  radial  stresses,  and  strains  are  generated 
when  the  stress  wave  travel  along  the  bar.  Such  variables  are  functions  of  radial  dimension,  so 
they  are  not  uniform  across  the  cross  section.  The  wave  dispersion  effect  exists,  which  means 
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that  the  harmonic  wave  with  higher  frequencies  travels  more  slowly  than  the  wave  with  lower 
frequencies.  Due  to  this  effect,  the  wave  shape  keeps  changing  as  it  travels  along  the  axial 
direction.  This  effect  should  be  considered  for  accurate  data  reduction. 

Many  authors  have  investigated  the  wave  dispersion  effect  in  order  to  modify  the  classical  1-D 
data  analysis  algorithm.  Their  study  is  based  on  the  work  of  Pochhammer  {4)  and  Chree  (5), 
who  independently  solved  the  problem  of  wave  propagation  in  an  infinite  cylindrical  rod  with 
finite  diameter  and  derived  the  frequency  equation.  This  frequency  equation  establishes  the 
relationship  between  phase  velocities  and  wavelengths  of  harmonic  waves.  Later  Davies  (2)  and 
Bancroft  (6)  solved  the  equation  numerically  for  some  limiting  cases.  With  the  help  of  modern 
computers,  the  frequency  equation  can  be  easily  solved  over  the  complete  range  of  geometric 
and  material  parameters.  Follansbee  and  Frantz  (7)  and  Gong  et  al.  (5)  investigated  the  effect  of 
wave  dispersion  correction  for  the  SHPB  test  and  found  that  significant  improvement  on  the 
stress-strain  data  can  be  achieved  by  wave  dispersion  correction.  During  dispersion  correction, 
Fourier  transformation  (FT)  is  used  to  decompose  the  measured  strain  data  into  its  harmonic 
waves.  The  dispersion  correction  methodology  is  then  applied  for  each  component  frequencies. 
Finally,  inverse  FT  is  performed  to  get  the  dispersed  wave  at  the  IB-S  and  S-TB  interfaces. 

Other  dispersion  correction  methods  were  also  suggested,  e.g.,  in  Li  and  Lambros  (P).  Currently 
wave  dispersion  correction  has  become  a  standard  procedure  in  SHPB  data  analysis. 

In  this  report,  the  wave  dispersion  in  cylindrical  tubes  for  the  SHPB  application  is  investigated. 
This  study  has  become  important  and  indispensable  because  recently  cylindrical  tubes  have  been 
used  for  SHPB  experiments.  For  example,  Dowling  et  al.  (10)  performed  a  dynamic  punch  shear 
test  to  study  the  mechanical  properties  of  metals.  This  technique  uses  a  transmission  tube  in 
place  of  the  TB,  such  that  the  IB  can  shear  a  cylindrical  disc  through  the  specimen.  This  test 
method  is  still  being  extensively  used  today.  Another  example  is  determining  the  mechanical 
properties  of  low-strength,  low-impedance  materials  under  dynamic  loading.  If  solid  TBs  are 
used  for  the  SHPB  test,  the  noise  level  is  so  high  that  proper  interpretation  of  the  transmission 
signal  becomes  impossible.  To  solve  this  problem,  some  researchers  have  used  visco-elastic  bars 
e.g.,  see  Zhao  et  al.  (11).  The  disadvantage  of  using  visco-elastic  bars  is  that  additional 
assumptions  have  to  be  made  about  the  visco-elastic  material  properties,  which  may  not 
represent  the  true  characteristics  of  the  bar  materials  at  high  strain  rates.  Furthermore,  the 
material  properties  are  functions  of  temperature,  moisture,  and  aging  conditions,  which  are 
difficult  to  control.  In  addition,  the  dispersion  correction  algorithm  now  becomes  very  extensive 
and  time  consuming  due  to  the  visco-elastic  nature  of  the  bar  materials.  Other  researchers  have 
put  forward  the  idea  of  using  a  hollow  cylindrical  tube  for  the  transmission  part  of  the  SHPB  test. 
Recently,  Chen  et  al.  (12)  used  IB  and  transmission  tube  made  of  high-strength  Al  alloy  to 
measure  the  dynamic  mechanical  properties  of  silicon  (Si)  rubber  materials.  In  this  way,  they 
obtained  higher  amplitude  transmission  signal  that  was  achievable  using  a  solid  bar.  The 
problem  is  that  they  did  not  perform  any  dispersion  correction  on  the  incident  and  transmission 
signals,  which  were  collected  from  midpoints  of  the  bar/tube.  They  used  pulse  shaper  to  filter 
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the  high-frequency  components  in  the  waveform  and  then  claimed  that  dispersion  correction  was 
not  necessary.  This  argument  is  not  correct  because  the  same  argument  can  be  applied  to  the 
traditional  SHPB  test,  which  means  if  pulse  shaper  is  used,  dispersion  correction  is  not  necessary 
any  longer.  It  is  not  clear  how  to  define  the  lower  limit  of  high-frequency  wave  components 
without  a  dispersion  correction  methodology.  The  use  of  a  pulse  shaper  can  improve  the  result 
but  cannot  eliminate  the  necessity  of  dispersion  correction.  It  is  still  necessary  to  investigate  the 
wave  propagation  in  a  cylindrical  tube  and  find  the  appropriate  dispersion  relations,  and  is  the 
objective  of  the  present  study. 

Wave  propagation  in  infinitely  long  tubes  has  been  discussed  by  several  authors,  as  early  as  in 
the  1920s.  Ghosh  {13)  studied  the  longitudinal  motions  of  a  hollow  cylinder  with  3-D  theory.  In 
1952,  Fay  {14)  investigated  this  problem  again  and  derived  the  frequency  equation.  But  there 
was  an  error  in  his  equation,  which  was  pointed  out  by  Herrmann  and  Mirsky  {15).  The  work  of 
Mirsky  and  Herrmann  (M-H  model)  is  considered  the  most  comprehensive  treatment  {15,  16). 
They  used  both  3-D  theory  and  shell  theory  to  study  the  wave  motion  in  cylindrical  tube.  Using 
3-D  theory,  they  derived  the  exact  frequency  equation,  and  their  shell  theory  considered 
membrane  effect,  bending  effect,  transverse  shear  deformation,  and  rotary  inertia.  They  claimed 
that  the  results  of  3-D  theory  and  shell  theory  were  very  close. 

This  study  is  based  on  the  work  of  Mirsky  and  Herrmann  {16).  First,  the  two  frequency 
equations  will  be  solved  again  to  make  sure  that  the  shell  theory  really  gives  a  solution  to  phase 
speed,  which  is  accurate  enough.  Then  the  results  of  M-H  shell  theory  will  be  applied  to  the 
SHPB  test.  It  will  be  shown  that  the  first  mode  solution  is  correct  for  the  SHPB  application.  A 
modified  dispersion  correction  algorithm  will  be  suggested  to  make  accurate  prediction  of  strain 
at  the  IB-S  and  at  the  specimen-transmission  tube  (S-TT)  interfaces.  Finally,  the  effects  of  the 
tube  geometry  on  the  accuracy  of  data  reduction  will  be  discussed. 


2.  Solution  of  M-H  Model 


As  mentioned  earlier,  Herrmann  and  Mirsky  {15,  16)  derived  the  exact  frequency  equation  for 
cylindrical  tubes  by  solving  the  3-D  problem  and  one  approximate  frequency  equation  through 
shell  theory.  The  3-D  theory  gave  the  exact  frequency  equation  in  cylindrical  coordinates 
(r,0,z)  {15): 

F30  (c,  h  /  A)  =  [Kjo  (p)Koj  (y)  +  (p)Kio  (y)  +  (8  /(jr^pyab)) 

+  FK„  (y)Koo  (P)  +  (1  /  F)K„  (P)Koo  (y)]  +  [(1  +  B)  V(Fy  ^ab)]K„  (p)K„  (y) 

-  [(1  +  B)  /(yab)]  [aKj ;  (y)Kio  (p)  +  hK,^  (y)Koi  (p)] 

-[(l  +  B)/(Fyab)][aK„(P)Kio(y)  +  bK„(p)Ko,(y)]  =  0,  (1) 
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where 


K„n  (z)  =  J„  (zb)Y„  (za)  -  J„  (za)Y^  (zb),  e,  =  ^{X  +  2p.)/pc^  =  ^[plp 


and 


1),B  =  ( 


1),F  = 


Py 


(2) 


In  equation  2,  a  is  the  wave  number  2n/A  and  Jn,,Yjj^  are  the  Bessel  funetions,  a  and  b  are  the 
inner  and  outer  radii  of  the  tube,  X  and  qare  Lame’s  eonstants,  p  is  tube  material  density, 
e  is  the  phase  speed  to  be  solved,  h  =  b  -  a  is  the  thiekness  of  the  tube  wall  and  A  is  the  wave 
length.  Mirsky  and  Herrmann  also  derived  one  approximate  frequeney  equation  based  on  shell 
theory  (76).  This  equation  is  very  lengthy,  so  it  is  not  listed  here.  Basieally,  this  equation  is  in 
following  format: 

Fsheii(s,5)  =  Ai(m)(s"  -M)"(s"  -k")"5’^  +[A2(m)s®  +  A3(m)sV  A4(m)s"  +A5(m)]5*’ 


+  [Ag(m)s''  +  A7(m)s^  +  Ag(m)]5''  -A9(m)[s^  -  Ajo(m)]5^  =  0 . 


(3) 


where  s  =  c/c3,  5  =  h/A,  Aj,,,  Ajq  and  M  are  constants  which  are  functions  of  tube  geometry 
parameter  m  =  h/R  (76),  v  is  the  material  Poisson’s  ratio,  k  is  the  shear  coefficient  in  M-H 
shell  model  which  is  related  to  Poisson’s  ratio  v  (75).  Here  R  =  (a  +  b)  / 2  is  the  mean  tube 
radius.  From  equation  3,  it  can  be  seen  that  phase  speed  is  a  function  of  both  m  and  5  . 

Equation  1  is  transcendental  in  nature  but  can  be  solved  with  the  help  of  mathematical  software 
such  as  MAPLE  or  MATHEMATICA.  Erom  equation  2,  it  can  be  seen  that  P  and/or  y  could 
become  imaginary.  When  c  <  c^  and  c  <  c^ ,  both  P  and  y  are  imaginary.  When  c^  <  c  <  c^ , 
y  is  real  but  P  is  imaginary.  Wherever  one  variable  becomes  imaginary,  Bessel  function 
Im ,  must  be  used  in  place  of  ,  Y^^  .  Such  replacement  guarantees  that  the  imaginary 

number  i  will  be  eliminated  from  the  equation  so  that  the  real  roots  can  be  obtained. 

Equation  3  can  be  solved  easily  using  the  Newton-Raphson  method.  There  are  four  real 
solutions.  When  the  wavelength  is  very  long,  5  approaches  zero  and  the  first  and  second 
solutions  of  s  approach  the  shear  coefficient  k  .  This  coefficient  k  leads  numerically  to  the 
same  phase  speed  as  the  limiting  solution  for  wave  speed  in  cylindrical  rods  when  the 
wavelength  approaches  infinity.  So  it  is  interesting  to  see  that  M-H  shell  model  explains  the 
physical  meaning  of  this  limiting  solution. 

The  solutions  of  equations  1,  2,  and  3  for  various  tube  geometries  are  shown  in  figure  2.  There 
are  infinite  solutions  for  the  frequency  equations.  Here,  only  the  first-mode  solutions  are  shown. 
The  following  material  constants  are  assumed:  E  =  206.09  GPa  ,  p  =  7850  kg  /m^ ,  and 
Poisson’s  ratio  v  =  0.30.  It  can  be  seen  that  the  two  models  provide  very  close  solutions.  The 
maximum  difference  occurs  at  the  “valley”  of  the  curves.  As  an  example,  for  m  =  0.250 , 3-D 
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Figure  2.  Solution  of  the  frequency  by  3-D  theory  and  shell  theory  for 
various  m  values. 


theory  gives  0/0^=  0.5720,  shell  theory  gives  0/0^=  0.5916,  and  the  error  is  ~3%.  This  error  is 
0.54%,  0.29%,  0.21%  for  m  =  0.548 ,  m  =  0.667  (2/3)  and  m  =  1.000 ,  respeetively.  The  error 
decreases  when  m  increases.  Figure  3  shows  the  solutions  by  shell  theory  for  various  m  values. 
When  m  =  0.033  (1/30),  the  phase  speed  decreases  very  quickly  over  the  range  of  low  §  .  With 
increasing  m,  the  phase  speed  drops  more  slowly  with  5  .  There  exists  one  minimum  (valley) 
point  for  all  the  curves,  after  which  speed  starts  to  increase.  It  is  obvious  that  all  the  curves 
approach  to  one  limiting  speed  when  5  goes  to  infinity.  For  the  material  properties  used  here, 
this  limiting  phase  speed  is  0.927  c^ . 

One  important  observation  from  phase  speed  solutions  is  that  the  general  claim  “wave 
components  with  lower  frequencies  travel  faster  than  wave  components  with  higher 
frequencies,”  which  is  always  true  for  the  cylindrical  rod,  does  not  always  hold  true  for  the 
cylindrical  tubes.  For  small  m  values,  the  phase  speed  first  drops  then  increases.  When  m 
becomes  higher,  the  claim  may  hold  for  the  whole  range  of  5  . 

On  the  other  hand,  in  this  study,  it  will  be  shown  that  the  phase  speed  will  not  fall  into  the  range 
where  phase  speed  increases  with  §  .  When  we  discuss  the  results  later,  for  simplicity,  we  still 
use  the  claim  that  the  lower  frequency  wave  moves  faster. 

Bancroft  (6)  solved  the  frequency  equation  for  the  cylindrical  rod.  He  gave  the  value  of  c,^  /c^, 
when  §  approaches  infinity.  Herrmann  and  Mirsky  (75)  also  gave  these  values.  The 
comparison  is  shown  in  figure  4.  It  is  clear  that  these  two  sets  of  solutions  are  exactly  the  same, 
which  also  verifies  the  correctness  of  M-H  models. 


5 


First  Mode 


Figure  3.  Solution  of  the  frequency  equation  by  shell  theory  for  various 
m  =h/R. 


Poisson's  Ratio,  v 

Figure  4.  C]^  /  C;.  when  5  =  h  /  A  approaches  infinity,  c^,  =  ^E/  p  . 

Now  it  is  proven  that  the  3-D  theory  and  the  shell  theory  give  solutions  that  can  be  considered 
the  same.  Shell  theory  solutions  will  be  used  for  further  investigation.  The  advantage  of  using 
the  shell  theory  model  is  that  the  frequency  equation  of  the  shell  model  is  not  transcendental,  so 
it  is  straightforward  to  include  the  model  into  a  computer  code. 

Mirsky  and  Herrmann  discussed  the  first-mode  solution  thoroughly.  But  they  were  not  sure  if 
wave  motion  in  tubes  always  follows  the  first-mode  solution,  and  they  did  not  mention  the 
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application  of  their  models  to  the  SHPB  experiment  either.  So  it  is  not  elear  if  the  first-mode 
solution  ean  be  used  to  determine  the  phase  speed  in  the  SHPB  experiments.  Follansbee  et  al. 
(7)  proved  that  the  solid  oylindrieal  bars  in  the  SHPB  experiment  vibrate  in  the  first  mode,  so  it 
may  also  be  true  for  the  hollow  oylindrieal  tubes.  To  prove  this  assumption,  3-D  finite  element 
analysis  is  performed  to  simulate  the  wave  motion  in  the  tubes  and  is  oompared  with  the  shell 
theory  for  the  first  mode. 


3.  Modification  of  Traditional  Dispersion  Algorithm 


3.1  Solution  of  Phase  Speed 

In  predioting  the  dispersed  waveform  by  M-H  shell  model,  the  first  step  is  to  determine  the  phase 
speed  of  every  harmonio  wave  oomponent.  Another  equation  is  required  in  addition  to  the 
e/e^  ~  h/ A  relationship  shown  in  figure  3: 

27r(e/eJ-(h/A)=  kcOo(h/eJ,  (4) 

whieh  is  equivalent  to 

o(k)  =  f  (k)  A  =  (kcOfl  /  27r)  •  A .  (5) 

Here,  k  =  l,,,K  =  N/  2  and  N  are  the  numbers  of  strain  data  oolleeted  by  the  strain  gages.  Phase 
speed  e(k)  ean  be  solved  by  finding  the  interseetion  of  the  two  eurves  represented  by  equations 
3-5  numerieally.  Feliee  (77)  performed  nonlinear  fitting  on  the  solution  of  frequeney  equation 
for  oylindrieal  rods  to  obtain  one  explioit  relationship  between  o/o^  and  h/A  .  Then  it  was  easy 
to  use  this  relationship  and  equation  4-5  together  to  determine  the  phase  speed.  For  oylindrieal 
tubes,  this  fitting  method  does  not  work  beoause  the  shape  of  the  o/o^  ~  h/  A  eurves  ohanges 
signifioantly  with  m  as  shown  in  figure  3.  On  the  other  hand,  the  nonlinear  fitting  is  not 
neoessary.  This  is  beoause  there  is  one  implioit  funotional  relationship  between  o/o^,  and  h/A 
in  equation  3.  Giving  one  h/A  value,  the  first-mode  solution  o/o^  ean  be  determined 
numerieally  from  equation  3.  With  the  help  of  this  observation,  the  phase  speed  of  wave 
oomponents  ean  be  determined  using  equations  3  and  4. 

3.2  Reconstruction  of  Waveform 

A  modification  is  suggested  in  this  study  on  dispersion  correction.  In  traditional  dispersion 
correction  algorithm,  the  following  equation  is  used  to  reconstruct  the  waveform; 

f  (t)  =  ^  +  X  [^k  cos(kcOot’ )  +  Bj^  sm(kcOot’ )] ,  (6) 

2  k=i 


with 
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t  =  t  -  Az 


^1  P 


v^k  ^oy 


(V) 


Here,  A^,  ,  and  Bj,are  the  coefficients  of  harmonic  wave  components  determined  by  Fourier 

transform  (FT).  Az  is  the  distance  between  the  strain  gage  and  the  destination  position.  Cj,  is 
the  phase  speed  of  the  wave  with  frequency  kcog .  This  modification  is  explained  in  figure  5. 
Suppose  there  are  two  sinusoidal  waves  sin(t)  and  sin(2t)  traveling  together  in  a  tube.  At 
location  A,  the  two  waves  are  in  phase.  Now  let  them  travel  forward  to  location  B.  Because  the 
wave  with  higher  frequency  travels  slower  than  the  wave  with  lower  frequency,  wave  sin(t)  will 
arrive  at  location  B  earlier  (at  time  T i)  than  sin(2t)  (at  time  T2).  From  time  T 1  to  time  T2,  the 
contribution  from  wave  sin(2t)  to  stress  and  strain  at  location  B  should  be  zero  because  it  has 
not  arrived  yet.  In  the  traditional  algorithm,  the  stress  and  strain  response  is  taken  to  be 
sin(2(t-  Lab  /c2+Lab/co)),  where  Lab  is  the  distance  between  location  A  and  location  B,  C2  is  the 
phase  speed  for  sin(2t),  and  Co  is  the  fundamental  speed. 


0  5  10  15  20 

Time 


Figure  5.  Explanation  of  new  dispersion  algorithm  for  SFIPB. 

Based  on  this  observation,  the  algorithm  for  inverse  fast  Fourier  transform  (FFT)  in  SHPB 
application  should  be  modified.  Suppose  the  initial  time  the  signal  arrives  at  one  certain  position 
is  t  =  0  .  At  time  t  =  tx ,  the  signal  should  be  reconstructed  as  in  the  following  pseudocode: 

Signal:  =Ao/2.0; 

For  k=l  to  N/2 

If  (tx<(  Az/c^-Az/Cq)) 
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Do  nothing; 


Else 

Signal <=  Signal+  A,^  cos( kcOg  (tx- zlz / c j,  +  Zlz / c ^ ))+  sin( kcOg  (tx-  Az! +  Az I Cq))', 

The  higher  frequeney  wave  travels  slower  than  the  lower  frequeney  wave.  The  lower  frequeney 
wave  reaches  location  B  at  time  Ti.  No  response  should  be  generated  by  the  higher  frequency 
wave  at  location  B  until  after  time  T2. 

After  the  coefficients  A^,  A^,  andBj,  are  determined,  this  modified  dispersion  algorithm  is 
applied  to  each  harmonic  component  in  reconstructing  the  waveform. 


4.  Finite  Element  (FE)  Analysis 


FE  analysis  was  performed  to  check  the  prediction  of  the  M-H  shell  model.  Stress  wave 
propagation  in  two  cylindrical  tubes  was  simulated.  These  tubes  have  m  =  h/R  values  of 
0.548-0.667.  The  dimensions  of  the  tube  with  m  =  0.548  are  as  follows:  inner  radius 
a  =  7.24  mm,  outer  radius  b  =  12.7  mm,  and  the  length  is  E  =  1.524  m.  Tubes  with  such 
dimensions  are  commercially  available  and  suitable  for  Hopkinson  bar  applications.  The 
dimensions  of  the  tube  with  m  =  0.667  are  as  follows:  inner  radius  a  =  5.46  mm,  outer  radius 
b  =  10.92  mm,  and  the  length  is  E  =  1.524  m.  Both  tubes  have  a  wall  thickness  of  5.46  mm. 

The  tube  material  is  modeled  as  steel  with  Young’s  modulus  E  =  206.09  GPa,  Poisson’s  ratio 
v  =  0.3 ,  and  density  p  =  7850  Kg/m  .  A  trapezoidal  stress  pulse  is  applied  at  one  end  of  the 
tube,  and  the  wave  propagation  in  the  tube  is  monitored.  The  other  end  is  assumed  to  be  traction 
free,  so  the  stress  wave  will  reflect  from  that  end.  Due  to  the  symmetric  nature  of  this  problem, 
only  one  quarter  of  the  tube  is  analyzed. 

Commercial  FE  analysis  code  ABAQUS/Explicit*  was  used  to  perform  the  analysis.  Element 
type  C3D8R,  which  is  8-node  brick  element  with  reduced  integration,  was  used  to  generate  the 
mesh.  The  numbers  of  elements  along  the  radial  direction,  tangential  direction,  and  axial 
direction  are  5,  20,  and  600,  respectively,  with  a  total  of  60,000  elements.  These  elements  are 
uniformly  distributed  along  the  axial  direction.  Figure  6  shows  the  model  and  the  mesh 
generation.  The  boundary  conditions  are  set  as  follows:  zero  x  displacement  along  edge  X  =  0 ; 
and  zero  y  displacement  along  edge  Y  =  0 .  Figure  7  defines  the  trapezoidal  stress  pulse  applied 
on  the  end  face  Z  =  0  .  The  stress  history  is  defined  with  several  time  and  stress  parameters  as 
described  in  figure  7.  The  values  used  in  the  present  analyses  are  as  follows:  t^  =0  ps, 
tg  =25  ps,  t(,  =35  ps,  tg  =  105  ps,  tg  =115  ps,  tp  =  140  ps  and  =500  MPa. 


ABAQUS/Explicit  is  a  registered  trademark  of  Abaqus  Inc. 
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Figure  7.  Definition  of  stress  wave  applied  on  the  end  surface  of  the 
cylindrical  tube. 


Figure  8  gives  the  strain  history  at  the  mid-length  of  tube  for  a  time  period  up  to  700  ps.  The 
results  of  FE  analyses  are  compared  with  M-H  shell  theory  and  are  presented  in  the  next  section. 
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Figure  8.  Strain  vs.  time  at  the  mid-length  of  the  tube. 


5.  Results  and  Discussion 


The  axial  strain  predictions  by  both  M-H  model  and  FE  analysis  for  the  two  cylindrical  tubes  at 
z  /  L  =  0.50-0.75  are  presented  in  figures  9  and  10.  When  the  M-H  model  is  used,  an  implicit 
assumption  is  made  that  the  wave  propagation  in  the  tubes  is  1-D,  and  strain  is  determined  from 
stress  by  Hooke’s  law  =  cr^z  /E  .  Plots  in  the  left  column  show  the  complete  strain  history, 
and  the  closer  views  in  the  peak  region  are  listed  in  the  right  column.  Erom  the  four  plots  in 
figures  9  and  10,  it  is  obvious  that  the  predictions  by  the  M-H  model  and  EEM  analysis  are  very 
close.  The  rising  parts  of  the  waves  are  almost  overlapping.  Eor  the  peak  part,  the  Pochhammer 
modes  match  very  well.  The  only  difference  is  that  for  the  M-H  model,  the  beginning  part  of  the 
waves  is  abnormal,  which  will  be  explained  later.  The  results  in  figures  9  and  10  prove  that  the 
cylindrical  tubes  in  the  SHPB  experiment  also  vibrate  in  the  first  mode. 

Eigure  1 1  compares  the  wave  dispersion  results  by  the  traditional  algorithm,  and  the  new 
algorithm  used  in  this  study.  Here  the  strain  data  are  plotted  in  the  opposite  sign.  The  results  by 
the  traditional  algorithm  are  not  smooth,  and  discontinuities  are  seen  over  the  peak  part  of  the 
waves.  On  the  contrary,  the  modified  algorithm  predicts  fairly  smooth  waveform.  Erom  the 
smooth  waves,  the  number  of  Pochhammer  modes  can  be  clearly  counted.  If  the  old  algorithm  is 
compared  to  PE  results,  it  is  nearly  impossible  to  reach  the  conclusion  that  wave  motion  in  tubes 
follows  the  first-mode  solution  of  M-H  model.  The  comparison  clearly  verifies  the  correctness 
and  necessity  of  the  new  algorithm. 
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Figure  9.  Axial  strain  predicted  by  M-Fl  model  and  FE  analysis  for  the  tube  with  m  =  0.548. 

But  why  does  the  traditional  algorithm  work  for  rods?  Figure  12  gives  an  example  for  the  use  of 
both  algorithms  on  the  rod  with  v  =  0.29.  It  is  seen  that  the  predicted  wave  by  traditional 
algorithm  still  has  some  discontinuous  response  in  the  later  part  of  the  wave  peak,  although  the 
difference  between  the  two  algorithms  is  not  as  significant  as  for  the  tubes.  From  figure  3,  it  can 
be  seen  that  the  phase  speed  in  rods  changes  very  slowly  with  6,  compared  with  phase  speed  in 
tubes.  This  leads  to  more  Pochhammer  modes  for  rods  than  the  tubes.  Because  the  irregularities 
occur  at  the  later  part  of  the  wave  peak,  the  traditional  algorithm  still  gives  wave  dispersion 
prediction  accurate  enough  for  engineering  application,  and  the  shortcoming  of  the  traditional 
algorithm  is  concealed.  In  this  study  of  wave  dispersion  in  tubes,  the  limitation  of  the  traditional 
wave  dispersion  algorithm  is  identified  first. 

The  first  part  of  the  wave  dispersion  prediction  by  the  new  algorithm  shows  abnormal  behavior. 
This  is  because  at  the  beginning,  relatively  few  wave  components  contribute  to  the  dispersed 
wave,  which  can  be  seen  from  the  algorithm  described  earlier.  The  M-H  model  assumes  infinite 
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(b)  z/L  =  0.50  ,  closer  view 
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(d)  z  /  L  =  0.75 ,  closer  view 


Figure  10.  Axial  strain  predicted  by  M-Fl  model  and  FE  analysis  for  the  tube  with  m  =  0.667. 


long  tubes  and  continuous  wave  components.  Under  these  assumptions,  at  any  instance,  all  the 
wave  components  contribute  to  the  whole  waveform,  so  this  abnormal  behavior  will  not  be  seen. 
This  phenomenon  cannot  be  seen  in  real  tests  or  in  an  FEM  simulation  either.  When  the  new 
algorithm  is  used  to  predict  the  dispersed  wave,  it  is  assumed  that  the  wave  has  period  of 
1400  |us  during  which  the  FT  is  performed.  This  input  wave  is  a  pulse  wave  rather  than  a 
continuous  wave.  The  new  algorithm  does  not  allow  any  summation  of  wave  components  if  the 
travel  time  condition  is  not  satisfied.  Thus  the  initial  data  are  spurious  till  all  the  wave 
components  have  arrived  at  the  particular  point  of  interest.  This  problem  can  be  resolved  by 
taking  a  wider  time  window  during  FT  and  reducing  the  window  during  reconstruction. 

When  rods  are  used  in  the  SHPB  test,  the  only  geometric  parameter  affecting  the  wave 
propagation  in  rods  is  the  radius.  If  tubes  are  used  in  the  SHPB  test,  two  more  parameters  are 
involved;  the  inner  radius  and  the  m  value.  Among  the  outer  radius,  inner  radius,  tube  wall 
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Figure  1 1 .  Comparison  of  new  and  traditional  dispersion  algorithms. 


Figure  12.  Dispersed  wave  in  rod  determined  by  the  new  algorithm  and 
traditional  algorithm.  Material  Poisson’s  ratio  =  0.29. 
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thickness  h  and  m,  only  two  are  independent.  The  diseussion  below  will  foeus  on  h  and  m.  It 
ean  be  seen  from  figure  3  that  when  m  takes  lower  values,  the  phase  speed  deereases  more 
quiekly  when  5  is  small.  For  a  fixed  m  value,  when  the  tube  wall  is  thinner,  the  phase  speed  for 
all  the  wave  eomponents  will  be  eloser  than  when  the  h  value  is  higher.  This  ean  be  seen 
graphieally  from  figure  3.  So  lower  m  values  and  higher  h  values  will  make  the  phase  speed 
deereases  faster  as  a  funetion  of  frequeney,  whieh  will  make  the  wave  beeome  wider  and  wider 
when  it  propagates  in  the  tube.  To  eateh  all  the  wave  eomponents,  a  wider  time  window  must  be 
used  to  proeess  the  strain  data.  Usually  a  eomputer  eode  is  used  to  automatieally  seleet  the 
proeessing  window.  For  example,  Leber’s  method  {18)  finds  the  proeessing  window  in  the 
following  way:  first  find  a  window  between  the  point  where  the  wave  rises  and  the  point  where 
the  wave  drops  to  zero  again.  Another  window  that  is  -33%  wider  than  the  first  window  and 
ineludes  the  first  window  in  the  middle  is  used  as  the  proeessing  window.  If  the  wave  beeomes 
too  wide,  an  ineorreet  window  may  be  ehosen  whieh  omits  some  important  eomponents. 
Furthermore,  the  ineident  wave,  and  the  refleeted  wave  will  not  eontain  the  same  amount  of 
information  beeause  the  refleeted  wave  is  wider  than  the  ineident  wave,  but  the  proeessing 
windows  are  still  the  same  width.  To  obtain  an  aeeurate  result,  it  is  reeommended  to  use  a  tube 
with  higher  m  value  and  lower  h  value.  Sueh  a  eombination  is  aetually  not  diffieult  to  realize. 

When  tubes  are  used  in  the  SHPB  teehnique,  several  eonfigurations  of  the  apparatus  are  possible. 
If  an  ineident  tube  and  a  transmission  tube  are  used  as  the  apparatus,  it  is  reeommended  that  one 
hollow  eylindrieal  speeimen  be  used.  This  requires  that  the  tubes  and  the  speeimen  have  the 
same  inner  and  outer  diameters.  The  advantage  of  this  eonfiguration  is  that  the  1-D  requirement 
of  the  SHPB  teehnique  ean  be  better  satisfied  than  using  rods.  Figure  13  shows  the  stress  history 
at  the  integration  point  of  all  the  five  elements  at  z  =  1 .00  m  along  the  radial  direetion.  It  is  seen 
that  for  all  five  elements,  the  axial  stress  history  vs.  time  is  nearly  the  same.  The  normal 
stresses  in  the  other  two  direetions  are  also  plotted.  Figure  13  also  shows  that  the  stress  a  is 
very  low  (maximum  is  13  MPa)  eompared  to  g^^  (maximum  is  580  MPa),  and  the  stress  g^^  is 
even  lower  (maximum  is  0.7  MPa).  So  the  stress  distribution  ean  be  eonsidered  as  1-D,  and  the 
axial  stress  distribution  is  uniform  aeross  the  eross-seetion  of  the  tube.  The  strains  in  x  and  y 
direetions  ean  be  eonsidered  only  due  to  Poisson’s  effeet.  Another  option  is  deseribed  by  Chen 
{12),  where  only  the  TT  is  used.  When  this  eonfiguration  is  used,  a  eap  must  be  inserted 
between  speeimen  and  TT  to  transfer  the  load.  It  is  not  elear  what  effeet  this  eap  might  have  on 
the  stress  wave  propagation  and  will  be  addressed  in  a  separate  work. 
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Figure  13.  Stress  history  of  all  elements  at  z  =  1.00  m  along  X  =  0  axis. 


6.  Conclusions 


In  this  report,  wave  dispersion  in  cylindrical  tubes  applied  to  the  SHPB  experiments  is  studied. 
The  main  conclusions  are  as  follows; 

1 .  The  M-H  shell  model  gives  the  same  result  as  the  3-D  model. 

2.  The  FE  study  shows  that  wave  propagation  in  tubes  follows  the  first-mode  solution  of  the 
M-H  shell  model. 

3.  A  new  dispersion  algorithm  is  suggested,  and  it  can  better  predict  the  dispersion  wave. 

4.  For  purpose  of  the  SHPB  application,  tubes  with  higher  m  and  lower  h  values  should  be 
used. 
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ENG  SYS  INC 
A  ALEXANDER 
13000  TENSOR  LANE  NE 
FLINTSTONE  MD  21530 

1  AAI CORP 

DR  N  B  MCNELLIS 
PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

1  OFC  DEPUTY  UNDER  SEC  DEFNS 
J  THOMPSON 

1745  JEFFERSON  DAVIS  HWY 
CRYSTAL  SQ  4  STL  501 
ARLINGTON  VA  22202 

3  ALLIANT  TECHSYSTEMS  INC 

J  CONDON 
ELYNAM 
J  GERHARD 
WVOl  16  STATE  RT  956 
PO  BOX  210 
ROCKET  CENTER  WV 
26726-0210 

1  PROJECTILE  TECHNOLOGY  INC 
515  GILES  ST 

HAVRE  DE  GRACE  MD  21078 

1  HEXCEL  INC 
RBOE 

PO  BOX  18748 

SALT  LAKE  CITY  UT  841 18 

1  PRATT  &  WHITNEY 
C  WATSON 

400  MAIN  ST  MS  114  37 
EAST  HARTFORD  CT  06108 


5  NORTHROP  GRUMMAN 
B  IRWIN 
K  EVANS 
D  EWART 

A  SHREKENHAMER 
J  MCGLYNN 
BLDG  160  DEPT  3700 
1 100  WEST  HOLLYVALE  ST 
AZUSA  CA  91701 

1  HERCULES  INC 
HERCULES  PLAZA 
WILMINGTON  DE  19894 

1  BRIGS  COMPANY 
J  BACKOFEN 
2668  PETERBOROUGH  ST 
HERNDON  VA  22071-2443 

1  ZERNOW  TECHNICAL  SERVICES 
L  ZERNOW 

425  W  BONITA  AVE  STE  208 
SAN  DIMAS  CA  91773 

1  GENERAL  DYNAMICS  OTS 
L  WHITMORE 

10101  NINTH  ST  NORTH 
ST  PETERSBURG  FL  33702 

2  GENERAL  DYNAMICS  OTS 
FLINCHBAUGH  DIV 

K  LINDE 
T  LYNCH 
PO  BOX  127 
RED  LION  PA  17356 

1  GKN  WESTLAND  AEROSPACE 
DOLDS 

450  MURDOCK  AVE 
MERIDEN  CT  06450-8324 

2  BOEING  ROTORCRAFT 
P  MINGURT 

P  HANDEL 

800  B  PUTNAM  BLVD 
WALLINGFORD  PA  19086 
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5  SIKORSKY  AIRCRAFT 
G  JACARUSO 
T  CARSTENSAN 
BKAY 

S  GARBO  MS  S330A 
J  ADELMANN 
6900  MAIN  ST 
PO  BOX  9729 

STRATFORD  CT  06497-9729 

1  AEROSPACE  CORP 
G  HAWKINS  M4  945 

2350  E  EL  SEGUNDO  BLVD 
EL  SEGUNDO  CA  90245 

2  CYTEC  FIBERITE 
MEIN 

W  WEB 

1440  N  KRAEMER  BLVD 
ANAHEIM  CA  92806 

2  UDLP 

G  THOMAS 
M  MACLEAN 
PO  BOX  58123 
SANTA  CLARA  CA  95052 

1  UDLP  WARREN  OFC 
ALEE 

3 1 20 1  CHICAGO  RD  SOUTH 
SUITE  B 1 02 
WARREN  MI  48093 

2  UDLP 

R  BRYNSVOLD 
P  JANKEMS  170 
4800  EAST  RIVER  RD 
MINNEAPOLIS  MN  55421-1498 

1  LOCKHEED  MARTIN 
SKUNK  WORKS 
D  FORTNEY 
ion  LOCKHEED  WAY 
PALMDALE  CA  93599-2502 

1  LOCKHEED  MARTIN 
R  FIELDS 
5537  PGA  BLVD 
SUITE  4516 
ORLANDO  FL  32839 


I  NORTHROP  GRUMMAN  CORP 
ELECTRONIC  SENSORS 
&  SYSTEMS  DIV 
E  SCHOCH  MS  V  16 
1 745 A  W  NURSERY  RD 
LINTHICUM  MD  21090 

1  GDLS  DIVISION 
D  BARTEL 

PO  BOX  1901 
WARREN  MI  48090 

2  GDLS 
DREES 
M  PASIK 
PO  BOX  2074 
WARREN  MI  48090-2074 

I  GDLS 

MUSKEGON  OPER 
M  SOIMAR 
76  GETTY  ST 
MUSKEGON  MI  49442 

I  GENERAL  DYNAMICS 
AMPHIBIOUS  SYS 
SURVIVABILITY  LEAD 
G  WALKER 
991  ANNAPOLIS  WAY 
WOODBRIDGE  VA  22I9I 

6  INST  FOR  ADVANCED 
TECH 
HFAIR 
I  MCNAB 
P  SULLIVAN 
S  BLESS 
W  REINECKE 
C  PERSAD 

3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 

I  ARROW  TECH  ASSOC 

1233  SHELBURNE  RD  STE  D8 
SOUTH  BURLINGTON  VT 
05403-7700 

I  R  EICHELBERGER 
CONSULTANT 
409  W  CATHERINE  ST 
BEL  AIR  MD2I0I4-36I3 
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SAIC 

G  CHRYSSOMALLIS 

8500  NORMANDALE  LAKE  BLVD 

SUITE  1610 

BLOOMINGTON  MN  55437-3828 

UCLA  MANE  DEPT  ENGR IV 
H  T  HAHN 

LOS  ANGELES  CA  90024-1597 

UNIV  OF  DAYTON 
RESEARCH  INST 
R  YKIM 
AKROY 

300  COLLEGE  PARK  AVE 
DAYTON  OH  45469-0168 

UMASS  LOWELL 
PLASTICS  DEPT 
N  SCHOTT 
1  UNIVERSITY  AVE 
LOWELL  MA  01854 

IIT  RESEARCH  CTR 
DROSE 
201  MILL  ST 
ROME  NY  13440-6916 

GA  TECH  RESEARCH  INST 
GA  INST  OF  TCHNLGY 
P  FRIEDERICH 
ATLANTA  GA  30392 

MICHIGAN  ST  UNIV 
MSM  DEPT 
R  AVERILL 
3515  EB 

EAST  LANSING  MI  48824-1226 

UNIV  OF  WYOMING 
D  ADAMS 
PO  BOX  3295 
LARAMIE  WY  82071 

PENN  STATE  UNIV 
R  S  ENGEL 

245  HAMMOND  BLDG 
UNIVERSITY  PARK  PA  1 680 1 


2  PENN  STATE  UNIV 
R  MCNITT 
C  BAKIS 

212  EARTH  ENGR 
SCIENCES  BLDG 
UNIVERSITY  PARK  PA  16802 

I  PURDUE  UNIV 

SCHOOL  OF  AERO  &  ASTRO 
CTSUN 

W  LAFAYETTE  IN  47907-1282 

I  STANFORD  UNIV 

DEPT  OF  AERONAUTICS 
&  AEROBALLISTICS 
STSAI 

DURANT  BLDG 
STANFORD  CA  94305 

I  UNIV  OF  MAINE 

ADV  STR  &  COMP  LAB 
R  LOPEZ  ANIDO 
5793  AEWC  BLDG 
ORONO  ME  04469-5793 

I  JOHNS  HOPKINS  UNIV 
APPLIED  PHYSICS  LAB 
P  WIENHOLD 
1 1 100  JOHNS  HOPKINS  RD 
LAUREL  MD  20723-6099 

I  UNIV  OF  DAYTON 
J  M  WHITNEY 
COLLEGE  PARK  AVE 
DAYTON  OH  45469-0240 

I  NORTH  CAROLINA  ST  UNIV 
CIVIL  ENGINEERING  DEPT 
W  RASDORF 
PO  BOX  7908 
RALEIGH  NC  27696-7908 

5  UNIV  OF  DELAWARE 

CTR  FOR  COMPOSITE  MTRLS 
J  GILLESPIE 
M  SANT ARE 
S  YARLAGADDA 
S  ADVANI 
D  HEIDER 
201  SPENCER  LAB 
NEWARK  DE  1 97 1 6 
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1  DEPT  OF  MTRLS 
SCIENCE  &  ENGRG 
UNIV  OF  ILLINOIS 
AT  URBANA  CHAMPAIGN 
J  ECONOMY 

1304  WEST  GREEN  ST  II5B 
URBANA  IE  6 1 80 1 

1  UNIV  OF  MARYLAND 

DEPT  OF  AEROSPACE  ENGRG 
A  J  VIZZINI 

COLLEGE  PARK  MD  20742 

1  DREXEL  UNIV 
A  S  D  WANG 
3141  CHESTNUT  ST 
PHILADELPHIA  PA  1 9 104 

3  UNIV  OF  TEXAS  AT  AUSTIN 
CTR  FOR  ELECTROMECHANICS 
J  PRICE 
A  WALLS 
J  KITZMILLER 
10100  BURNET  RD 
AUSTIN  TX  78758-4497 

3  VA  POLYTECHNICAL 

INST  &  STATE  UNIV 
DEPT  OF  ESM 
M  W  HYER 
K  REIFSNIDER 
R  JONES 

BLACKSBURG  VA  24061-0219 

1  SOUTHWEST  RESEARCH  INST 
ENGR  &  MATE  SCIENCES  DIV 
J  RIEGEL 

6220  CULEBRA  RD 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78228-0510 

1  BATELLE  NATICK  OPERS 
B  HALPIN 
313  SPEENST 
NATICK  MA  01760 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARE  WM  MB 
A  FRYDMAN 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 


ABERDEEN  PROVING  GROUND 

I  US  ARMY  ATC 
CSTE  DTC  AT  AC  I 
W  C  FRAZER 
400  COLLERAN  RD 
APG  MD  21005-5059 

91  DIRUSARL 

AMSRD  ARE  Cl 
AMSRD  ARE  O  AP  EG 
M  ADAMSON 
AMSRD  ARE  SL  BA 
AMSRD  ARE  SL  BB 
DBELY 

AMSRD  ARE  WM 
J  SMITH 
H  WALLACE 
AMSRD  ARE  WM  B 
A  HORST 
T  KOGLER 
AMSRD  ARE  WM  BA 
DLYON 

AMSRD  ARE  WM  BC 
J  NEWILL 
P  PLOSTINS 
A  ZIELINSKI 
AMSRD  ARE  WM  BD 
P  CONROY 
B  FORCH 
M  LEADORE 
C  LEVERITT 
RLIEB 

R  PESCE  RODRIGUEZ 
BRICE 

AMSRD  ARE  WM  BF 
S  WILKERSON 
AMSRD  ARE  WM  M 
BFINK 
J  MCCAULEY 
AMSRD  ARE  WM  MA 
L  GHIORSE 
S  MCKNIGHT 
E  WETZEL 
AMSRD  ARE  WM  MB 
J  BENDER 
T  BOGETTI 
L  BURTON 
R  CARTER 
KCHO 

W  DE  ROSSET 
G  DEWING 
R  DOWDING 
W  DRYSDALE 
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R  EMERSON 

AMSRD  ARE  WM  TB 

D  HENRY 

P  BAKER 

D  HOPKINS 

AMSRD  ARE  WM  TC 

RKASTE 

R  COATES 

L  KECSKES 

AMSRD  ARE  WM  TD 

M  MINNICINO 

D  DANDEKAR 

B  POWERS 

T  HADUCH 

D  SNOHA 

T  MOYNIHAN 

J  SOUTH 

M  RAFTENBERG 

M  STAKER 

S  SCHOENFELD 

JSWAB 

T  WEERASOORIYA 

J  TZENG 

AMSRD  ARE  WM  TE 

AMSRD  ARE  WM  MC 

A  NIILER 

J  BEATTY 

J  POWELL 

R  BOSSOLI 
ECHIN 

S  CORNELISON 
D  GRANVILLE 
BHART 
J  LASALVIA 
J  MONTGOMERY 
F  PIERCE 
ERIGAS 
W  SPURGEON 

AMSRD  ARE  WM  MD 
B  CHEESEMAN 
P  DEHMER 
R  DOOLEY 
G  GAZONAS 
S  GHIORSE 
C  HOPPLE 
M  KLUSEWITZ 
WROY 
J  SANDS 
D  SPAGNUOLO 
S  WALSH 
S  WOLF 

AMSRD  ARE  WM  RP 
J  BORNSTEIN 
C  SHOEMAKER 

AMSRD  ARE  WM  T 
B  BURNS 

AMSRD  ARE  WM  TA 
W  BRUCHEY 
M  BURKINS 
W  GILLICH 
B  GOOCH 
T  HAVEL 
E  HORWATH 
M  NORMANDIA 
J  RUNYEON 
M  ZOLTOSKI 
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1  LTD 

R  MARTIN 
MERE 

TAMWORTH  RD 
HERTFORD  SGI 3  7DG 
UK 

1  SMC  SCOTLAND 
P  WLAY 
DERA  ROSYTH 
ROSYTH  ROYAL  DOCKYARD 
DUNFERMLINE  FIFE  KY  1 1  2XR 
UK 

1  CIVIL  AVIATION 
ADMINSTRATION 
T  GOTTESMAN 
PO  BOX  8 

BEN  GURION  INTRNL  AIRPORT 

LOD  70150 

ISRAEL 

1  AEROSPATIALE 
S ANDRE 

ABTE  CCRTE  MDI32 
316  ROUTE  DE  BAYONNE 
TOULOUSE  31060 
FRANCE 

1  DRA  FORT  HALSTEAD 
P  N  JONES 

SEVEN  OAKS  KENT  TN  I47BP 
UK 

1  SWISS  FEDERAL  ARMAMENTS 
WKS 
WLANZ 

ALLMENDSTRASSE  86 
3602  THUN 
SWITZERLAND 

1  DYNAMEC  RESEARCH  LAB 
AKE  PERSSON 
BOX  201 

SE  I5I  23  SODERTALJE 
SWEDEN 


I  ISRAEL  INST  OF  TECHLGY 
S  BODNER 

FACULTY  OF  MECHANICAL 

ENGR 

HAIFA  3200 

ISRAEL 

I  DSTO 

WEAPONS  SYSTEMS  DIVISION 
N  BURMAN  RLLWS 
SALISBURY 

SOUTH  AUSTRALIA  5108 
AUSTRALIA 

I  DBF  RES  ESTABLISHMENT 
VALCARTIER 
A  DUPUIS 

2459  BLVD  PIE  XI  NORTH 
VALCARTIER  QUEBEC 
CANADA 

PO  BOX  8800  COURCELETTE 
GOA  IRO  QUEBEC 
CANADA 

I  ECOLE  POLYTECH 
J  MANSON 
DMX  ETC 

CH  1015  LAUSANNE 
SWITZERLAND 

1  TNO  DEFENSE  RESEARCH 
R  IJSSELSTEIN 
ACCOUNT  DIRECTOR 
R&D  ARMEE 

PO  BOX  6006 
2600  JA  DELFT 
THE  NETHERLANDS 

2  FOA  NATL  DEFENSE  RESEARCH 
ESTAB 

DIR  DEPT  OF  WEAPONS  & 
PROTECTION 
B  JANZON 
R  HOLMLIN 
S  172  90  STOCKHOLM 
SWEDEN 
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2  DEFENSE  TECH  &  PROC 
AGENCY  GROUND 
I  CREWTHER 
GENERAL  HERZOG  HAUS 
3602  THUN 
SWITZERLAND 

1  MINISTRY  OF  DEFENCE 
RAFAEL 

ARMAMENT  DEVELOPMENT 
AUTH 

M  MAYSELESS 
PO  BOX  2250 
HAIFA  3 1 02 1 
ISRAEL 

1  TNO  DEFENSE  RESEARCH 
I H  PASMAN 
POSTBUS  6006 
2600  JA  DELFT 
THE  NETHERLANDS 

1  B  HIRSCH 

TACHKEMONY  ST  6 

NETAMUA426II 

ISRAEL 

1  DEUTSCHE  AEROSPACE  AG 
DYNAMICS  SYSTEMS 
MHELD 
PO  BOX  1340 

D  86523  SCHROBENHAUSEN 
GERMANY 
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